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We apply the energy-momentum tensor which is coordinate independent to calculate the 
energy content of the axisymmetric solutions. Our results are compared with what have 
been obtained before within the framework of Einstein general relativity and Mpller’s tetrad 
theory of gravitation. 


1. Introduction 


General geometric arena of PGT, the Riemann-Gartan space f/4, may be a priori restricted 
by imposing certain conditions on the curvature and the torsion. Thus, Einstein’s GR is 
defined in Riemann space V4, which is obtained from f/4 by the requirement of vanishing 
torsion. Another interesting limit of PGT is the teleparallel or Weitzenbock geometry T4. 
The vanishing of the curvature means that parallel transport is path independent. The 
teleparallel geometry is, in sense, complementary to Riemannian: curvature vanishes, and 
torsion remains to characterize the parallel transport. Of particular importance for the 
physical interpretation of the teleparallel geometry the fact that there is a one-parameter 
family of teleparallel Lagrangians which is empirically equivalent to GR [8, 23, 9, 12]. For the 
parameter value B = 1/2 the Lagrangian of the theory coincides, modulo a four-divergence, 
with the Einstein-Hilbert Lagrangian, and defines (TEGR). 

The teleparallel equivalent of general relativity (TEGR) is a viable alternative geomet¬ 
rical description of Einstein’s general relativity written in terms of the tetrad field [10] and 
continue to be object of thorough investigations [33, 34, 35, 36]. In the framework of the 
TEGR it has been possible to address the longstanding problem of defining the energy, mo¬ 
mentum and angular momentum of the gravitational field [13, 14, 15]. The tetrad field seems 
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to be a suitable field quantity to address this problem, because it yields the gravitational 
field and at the same time establishes a class of reference frames in space-time [16]. Moreover 
there are simple and clear indications that the gravitational energy-momentum defined in 
the context of the TEGR provides a unified picture of the concept of mass-energy in special 
and general relativity. 

The tetrad formulation of gravitation was considered by Mpller in connection with at¬ 
tempts to define the energy of gravitational field [5, 6]. For a satisfactory description of 
the total energy of an isolated system it is necessary that the energy-density of the gravita¬ 
tional field is given in terms of first- and/or second-order derivatives of the gravitational field 
variables. It is well-known that there exists no covariant, nontrivial expression constructed 
out of the metric tensor. However, covariant expressions that contain a quadratic form of 
first-order derivatives of the tetrad field are feasible. Thus it is legitimate to conjecture that 
the difficulties regarding the problem of defining the gravitational energy-momentum are 
related to the geometrical description of the gravitational field rather than are an intrinsic 
drawback of the theory [13, 17]. 

Definition of the angular momentum of the gravitational field is given in the framework 
of the TEGR [18]. In similarity to the definition of the gravitational energy-momentum, 
Maluf et al. [18] have interpreted the appropriate constraint equations as equations that 
defined the gravitational angular momentum. This definition turned out to be coordinate 
independent. The definition of is invariant under global SO{3, 1) transformations. It has 
been argued elsewhere [19] that it makes sense to have a dependence of on the frame. 
The energy-momentum in classical theories of particles and fields does depend on the frame, 
and it has been asserted that such dependence is a natural property of the gravitational 
energy-momentum. The total energy of a relativistic body, depends on the frame. 

A well posed and mathematically consistence expression for the gravitational energy has 
been developed [13]. It arises in the realm of the Hamiltonian formulation of the TEGR [20] 
and satisfies several crucial requirements for any acceptable definition of gravitational energy. 
The gravitational energy-momentum P“ [13, 14] obtained in the framework of the TEGR 
has been investigated in the context of several distinct configuration of the gravitational 
filed. For asymptotically fiat space-times P*-°^ yields the ADM energy [31]. In the context 
of tetrad theories of gravity, asymptotically fiat space-times may be characterized by the 
asymptotic boundary condition 

1 

Ca/i = T 

and by the condition = 0(l/r^) in the asymptotic limit r —)■ cx), with 

rjab = (—1,-|-1,-|-1,-|-1) is the metric of Minkowski space-time. An important property of 
tetrad fields that satisfy the above as that in the fiat space-time limit one has x, y, z) = 
6°'^, and therefore the torsion tensor = 0. Maluf [16] has extended the definition of 
P“ for the gravitational energy-momentum [13, 20] to any arbitrary tetrad fields, i.e., for 
the tetrad fields that satisfy 7^ 0 for the fiat space-time. The redefinition is the only 
possible consistent extension of P“, valid for the tetrad fields that do not satisfy the above 
equation. 

Recently, Sharif and Amir (2007) have found the teleparallel (TP) version of the non¬ 
null Einstein Maxwell solutions [27]. Then, they have used the TP version of Mpller (1978) 
to evaluate the energy-momentum distribution of these solutions. They have found that 
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the energy content in the TP theory is equal to the energy in GR (as found by Sharif and 
Fatima (2006)) plus some additional part. Also they have discussed three possibilities for 
the axial-vector held. 

It is the aim of the present work to calculate the energy content of the axisymmetric 
solutions using the dehnition of the gravitational energy which is coordinate independent. 
In Sect. 2, a brief review of the derivation of the held equations of the gravitational held 
is given. A summary of the derivation of energy and angular momentum in TEGR is also 
given in Sect. 2. In Sect. 3, we derive the axially symmetric solutions in TEGR and then, 
calculate their energy content. The hnal section is devoted to discussion and conclusion. 


2. The TEGR for gravitation 


In a spacetime with absolute parallelism the parallel vector held dehne the nonsym- 
metric affine connection 


■pA drf. 

/if/ 


e/e“ 




( 1 ) 


where = duCap*. The curvature tensor dehned by T^^y is identically vanishing, however. 
The metric tensor is given by 

9^u fif' ui ( 2 ) 


with the Minkowski metric rjab = diag(-|-l , —1 , —1 , —1) h 


The Lagrangian density for the gravitational held in the TEGR, in the presence of matter 
helds, is given by^ [13, 16] 


Cg = eLc 


167r V 4 2 



e 

167r 


S abcrp 

abc 


( 3 ) 


where e = det{e°‘^). The tensor is dehned by 


'Y^ahc <^6f. 





and T“ are the torsion tensor and the basic vector held dehned by 


( 4 ) 


T“ def-ga t^a = Q - d 

flU ^ X-‘- flU U 111 

and 

def. rpv n rpa def. 

*spacetime indices g, v, - ■ ■ and SO(3,l) indices a, b • • • run from 0 to 3. Time and space indices are 
indicated to /x = 0 , i, and a= ( 0 ), (i). 

t Latin indices are rasing and lowering with the aid of gab and 77 “**. 
throughout this paper we use the relativistic units , c = G = 1 and k = Stt. 


( 5 ) 

( 6 ) 
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The quadratic combination T,°-’^'^Tahc is proportional to the scalar curvature -R(e), except for 
a total divergence term [17]. represents the Lagrangian density for matter helds. 

The gravitational held equations for the system described by Lq are the following 


where 




= eT 

— afi' 


It is possible to prove by explicit calculations that the left hand side of the symmetric part 
of the held equations (7) is exactly given by [13] 


e 
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The axial-vector part of the torsion tensor is dehned by 


A -e 


jumper upa 

2 ^pvpcr I ) 


where 


V Q^pupcT’) 


( 8 ) 


with 7 ,^po- = rj^-^eauebp ; a being the contorsion tensor and 6^upa is completely antisymmetric 
and normalized as (5oi23 = — 1- 

In the context of Einstein’s general relativity, rotational phenomena is certainly not a 
completely understood issue. The prominent manifestation of a purely relativistic rotation 
ehect is the dragging of inertial frames. If the angular momentum of the gravitational held 
of isolated system has a meaningful notion, then it is reasonable to expect the latter to be 
somehow related to the rotational motion of the physical sources. 

The angular momentum of the gravitational held has been addressed in the literature 
by means of diherent approaches. The oldest approach is based on pseudotensors [29, 26], 
out of which angular momentum superpotentials are constructed. An alternative approach 
assumes the existence of certain Killing vector helds that allow the construction of conserved 
integral quantities [1]. Finally, the gravitational angular momentum can also be considered 
in the context of Poincare gauge theories of gravity [24], either in the Lagrangian or in the 
Hamiltonian formulation. In the latter case it is required that the generators of spatial rota¬ 
tions at inhnity have a well dehned functional derivatives. From this requirement a certain 
surface integral arises, whose value is interpreted as the gravitational angular momentum. 

The Hamiltonian formulation of TEGR is obtained by establishing the phase space vari¬ 
ables. The Lagrangian density does not contain the time derivative of the tetrad component 
Cao- Therefore, this quantity will arise as a Lagrange multiplier [30]. The momentum canon¬ 
ically conjugated to Cq* is given by = SL/Seai- The Hamiltonian formulation is obtained 
by rewriting the Lagrangian density in the form L = pq — H, in terms of eai,n“* and the 
Lagrange multipliers. The Legendre transformation can be successfully carried out and the 
hnal form of the Hamiltonian density has the form [20] 


i/ = e„oC'“ + a,fc^'= + /3fc^^ 


(9) 
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plus a surface term. Here aik and I3k are Lagrange multipliers that are identified as 


c^ik 2 and 




OOk, 


( 10 ) 


and (7“, r*^ and L^ are first class constraints. The Poisson brackets between any two held 
quantities F and G is given by 


{F, G} = I d^x 


6F 6G 
6eai{x) (5n“*(a;) 


6F 6G \ 
5n“*(a;) Scaiix)) ' 


( 11 ) 


We recall that the Poisson brackets |P*-^ (a;), P^^(a;)| reproduce the angular momentum alge¬ 
bra [17]. 

The constraint G°' is written as G^ = + /i“, where is an intricate expression 

of the held variables. The integral form of the constraint equation C" = 0 motivates the 
dehnition of the gravitational energy-momentum four-vector [17] 

= - / d^xdiW\ ( 12 ) 

Jv 

where V is an arbitrary volume of the three-dimensional space. In the conhguration space 
we have 

(13) 

K 

with 

= + where ^ . 

The emergence of total divergences in the form of scalar or vector densities is possible in 
the framework of theories constructed out of the torsion tensor. Metric theories of gravity 
do not share this feature. By making A = 0 in Eq. (13) and identifying H"* in the left side 
of the latter, the integral form of Eq. (13) is written as 

P“ = . (14) 

Eq. (14) suggests that P“ is now understood as the gravitational energy-momentum [19]. 
The spatial component P^*^ form a total three-momentum, while temporal component 
is the total energy [26]. 

It is possible to rewrite the Hamiltonian density of Eq. (9) in the equivalent form [18] 

H = eaoG^ + -AafeP^^, with \ab = —Aba, (15) 
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are the Lagrangian multipliers that are identified as = aik and Aqa; = — A^o = Pk- The 
constraints [20] embodies both constraints and F^ by means of the relation 

= e„*e^,^F“^ and F^ = F°^ = ea°e^,^F“^ (16) 

The constraint F“^ can be reads as 

pafe _ ^ab ^ . (17) 

K ^ ' 

In similarity to the definition of P“, the integral form of the constraint equation F“^ = 0 
motivates the new definition of the space-time angular momentum. The equation F“^ = 0 
implies 

. (18) 

K ^ ^ 

Maluf et ah [17, 18] defined 

= 2 / (19) 

JV 

as the 4-angular-momentum of the gravitational field for an arbitrary volume V of the three- 
dimensional space. In Einstein-Cartan type theories there also appear constraints that sat¬ 
isfy the Poisson bracket given by Eq. (11). However, such constraints arise in the form 
n[*d = 0, and so a definition similar to Eq. (19), i.e., interpreting the constraint equation as 
an equation for the angular momentum of the field, is not possible. Definition (19) is three- 
dimensional integral. The quantities P“ and are separately invariant under general co¬ 
ordinate transformations of the three-dimensional space and under time reparametrizations, 
which is an expected feature since these definitions arise in the Hamiltonian formulation of 
the theory. Moreover, these quantities transform covariantly under global SO{3, 1) transfor¬ 
mations [18]. 


3. Energy content of axisymmetric solutions 


Now we are going to calculate the energy content of the axisymmetric tetrad field that 
has the form [27] 


I 


fe/) = 


B{p,z) . 

B{p ,z) 


0 


o-K(p,z) 


0 


cos ( 


F(p) 

0 


cos^ 


0 


np) 

1 


sin0 


0 \ 
0 


np) 


coscj) 0 


( 20 ) 


0 / 
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Using Eq. (2) the associated metric of the tetrad held given by Eq. (20) takes the well 
known form 

ds"^ = dt^ — — B‘^{p , z)^ dcj)^ — + 2B{p , z)dtd((), (21) 

B{p ,z), K{p, z) and F{p) are unknown functions which satisfy the following relations 
B = FW\ B' = ~aF (fU^ _ ^ 

K'= ~aFWW\ fU + FF-ifU + fU" = 0, (22) 

where dot and prime denoting the derivatives w.r.t. p and respectively. Here a is a constant 
and fU is an arbitrary function of p and 2 : in general [27]. McIntosh’s give a solution in the 
form W = —2bz while McLenaghan et. ah solution’s has the form W = 2lnp [2]. The above 
metric represents a hve classes of non-null electromagnetic held and prefect huid solutions 
which possesses a metric symmetry not inherited by the electromagnetic held and admits a 
homothetic vector held. Two of these classes contain electrovac solutions as special cases, 
while the other three necessarily contain some huid. Generalization of metric given by Eq. 

( 21 ) is given in [ 11 ]. 

Applying the tetrad held of Eq. (20) to the held equations (7) we get the non-vanishing 
components to have the form 

_„-2K 

T% = ^3 (2B[FB + FB" - BF] + ?>F[B‘^ + B'^] - AF‘^[Fk + FK” + F]) , 

_p-2K 

T °2 = - r— {f\FB + FB” -BF- 2BF] - B^[BF - FB - FB”] + 2BF[B^ + B'^]) , 

327rr'^ ^ ^ 

-2K -2K 

-2K -2K 

T% = ^3 (f[F + B”] - BF) , T^2 = ^3 (f[B‘^ + 2BB + B’^ + 2BB”] - 2BBF + AF^[k + K”]) , 

-2K 

T\ = {b'^ -B^ + AF[F - Fk]) , (23) 

where is the energy-momentum tensor. 

A special solution of the above non linear P.D.E can be obtained by choosing [27] 

777 

B{p, z) = -e””, F{p) = e””, K{p, z) = 0, (24) 

n 

where m and n are constants. The above solution reproduce the well known solution which 
is known as the electromagnetic generalization of the Gddel solution [27, 2 ]. 

The second solution can be obtained by choosing 

pCip 

B(p,z)=e‘“-, = A'(p.z) = 0, (25) 

which known as the space time homogenous Gddel metric [27, 2 ]. 
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Now let us calculate the energy content of the tetrad (20) using (12). To do so let us 
calculate the non-vanishing components of the torsion tensor. Using Eq. (20) in Eq. (5) we 
get 

= j; ((e^-F)[e^BFsin^</) +B(Fcos </) + !)]-FB) , = -e^^BK'sin </)cos </) 

T (°^23 = B', = sin0(F-e^), T«i 3 = cos0, = (e^ - B)(e^ sin^ (/> + cos (/>), 

T( 2)^3 = -e^^K'sine/) cos (/), = -e^K. (26) 

The non-vanishing components of the basic vector are 

Ti = (F-e^)(e^sin2 0 + cos0)+e^h:, T 2 = sin0(F-e^), T 3 = K'coscj). (27) 


Pereira et al. [21] have proved that the axial vector tensor plays the role of the gravit- 
omagnetic component of the gravitational field in the case of slow rotation and weak field 
approximations. The non-vanishing components of the axial vector tensor, dehned by Eq. 
(8), associated with the tetrad held given by Eq. (20) are 




K' sin 0 cos 0(E^ — 2B‘^) 


3F 


A^ = 


^-2K 

3F 


A^ = 


2BK' sin 0 cos ( 


3F2 


(^2BF{l — e ^F){sin^4) + e ^cos0}-|-5e ^ 


3F 
+ BF)] . 


(28) 


It is of interest to compare our results with that obtained before by Sharif and Amir (2007). 
They have calculated the axial vector of the tetrad given by Eq. (20) and wrote the non¬ 
vanishing components as A^ which coincides with what we have obtained in Eq. (28) and the 
other component as A^ that has the form 


3F 

which is completely different from that obtained in Eq. (28). Therefore, the analysis related 
to the axial vector part given in Ref. [27] will completely now need some modifications which 
we will do. The spacelike axial vector can now be written [27] 



A — ^/—giiA^ep -I- ^/--gffA?!^ + \/—g^A^ez, 


(30) 


where ip, iff, and iz are unit vectors along the radial p, 0 and z-directions respectively. Using 
Eq. (28) in (30) we get 




A = 


B' ^ 2^/F^ - B^BK' 


3F 


C/j + 


sm 0 cos ( 


3F 


"60 


+ 


3F2 


1 —e ^F) 2i?F{sin^ 0 + e cos 0}-|-Fe 


.,-K 


e-^^FB]ez. (31) 


The spin precession of a Dirac particles in teleparallel gravity is related to the axial vector 
by [27, 25] 


dS 

dt 


-b X S, 


(32) 



where S is the spin vector of a Dirac particles and b = 3/2A where A is given by Eq. (31). 

The direct evidence for the coupling of intrinsic spin to the rotation of the Earth has 
become available [3]. According to the TEGR every spin | particle in the laboratory has 
an additional interaction Hamiltonian. However, such intrinsic spin must precess in a sense 
opposite to the sense of rotation of the Earth as measured by the observer. The corresponding 
extra Hamiltonian associated with such motion would be of the form [4] 

6H = -h- a, (33) 


where a is the spin of the particle. 

To calculate the energy density associated with the tetrad field given by Eq. (20) we 
must calculate which dehned in Eq. (4). The necessary non-vanishing components of 


{1 - e-^F){F{sm^ (p + e-^ cos(j)){B^ - F^) + e-^- e-'^^BFB + e-^F^kj , 

-e-^{e-^BB' - cos cpF^K') 


^002 _ 


2F3 

\ 

— sin (pk^ — F) 


^003 _ 


2F2 


2F2 




^101 ^ -Hsin 0(e “ - e 




1)F) 


2F2 


^102 _ 


4F3 


1 - e-^F){2BF[sik 0 + cos 0] + e-^B} - e-^^FB) , 

-K 


^103 _ 1 -2if^^/sin0cos0, = „-K 

4 AF'^ 


^202 _ Bsm<p 


- E)) , = 


^302 _ 


—e 


2F2 
-2Kb' 


203 


j-2iC 

AF‘^ ’ 


{B - Be-^F -e-^FB), 

e~^^BK' sin 6 cos ( 


J.301 _ 


- F). 


AF‘^ ’ 2F2 

Using (34) in (13) we get 


(34) 




d 


—(^1 

kF‘^ \ [dp 

- cos (pF^^iK'e^) + 
dz 


{BB) + ^{BB') + sik(P 
dz 


B^ -F^ 


d_ 

dp 


(Fe^) \ + BBF [1 + 2F cos (p] + B^ 


FF - 2F^ 


! 

dp 


2K\ K: 


cos0^(FF) 

dp ' 

1 -I- F cos^ 0| 

{f^ + f^} 

F2 + F^f] - F^ cos 0 

-^{Fc^ 1-1- cos 0e^ } 
dp '' L J J 


dp \B 


d 


+ 2 cos (pk-{e^B'^F) 
dp 


-F^^{ke^) Fe^Bsiki 
dp 


2FFB - {fBK + 2FB - f} 


-h e^^F cos I 


1 - e-^F 


( 35 ) 


The non-vanishing components of the momentum density have the form 

^( 1)0 ^ sm0cos0e /d_ ^ Fe^^F + FF^ ^ (e^^k) - Be^^F^ 

kF-^ \ dp ^ dp dp ^ 
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-F^{BF) + Fe'^^BF - 2e^^FB cos0 + 2FBe^ cos (j)F + FBe^ - FBF + — cos (j)-^{BK') 
dp 2 dz 


n(2)o 

]^(3)0 


COS 0 


kF2 

\ 

2 sin 6 


Be^F - BF^ - ^{Be^k - B + BF + FB + FB”) + F‘^Be^^ cos0 - F'^Be^^ cos0F 


kF 

\ 

-Be^^{B' + AK' 


Fe^^[B' + ?,BK'] - ^^^^^ {FBKK' + BFK' + FBK' + FBK') 


(36) 


It is of interest to note that if (-F(p) = const, K{p,z) = Ki{p), B{p,z) = B{p)) then the 
component of the momentum density = 0. 

Now let us repeat the same calculations using the solution given by Eq.(24), in this case 
the basic vector has the non-vanishing components 


Ti = {ne"''^ — 1) (sin^ 0 -1- cos 0), T 2 = sin (f){ne^^ ~ 1), 


(37) 


and the non-vanishing components of the axial vector held are 


= 


A= - 
3 


m 


m 


2-^(1 — ne”^){sim 0 -I- COS0} -I-e 

n n 


777 777 

2—(1 — ne'^^){siik 0 -|- cos 0} -I- 


— 2m 


n n 

The corresponding extra Hamiltonian [7] is given by 


2m 1 Cz- 


(38) 


(39) 


SH = —b ■ a = -- 


m , 


m 


2—(1 — ne'^^){sin^ 0-1-COS0}-I-e 

n n 


— 2m ) 62 ■ a, 


and the component of the energy density is given by 


„np 


(40) 


n(o)o = — [n^ — rk) [sin^ 0 -|- cos 0] (1 — 2ne”^) -|- ne cos 0{e — n} + 2m^n .(41) 


UK 


The non-vanishing components of the momentum density have the form 
n'o» = (,-.{1 _ 3 ^ _ „} + £ T!(i _ + „(3 cos^ ^ -1)' 

Thhx/ \ ^ 


n(2)0 _ - ^ ^2 f^ 3^g2np _ ^ _ ^^np ^ ^e^p ^^^2 ^ ^ 3 ^ 0e2^^(l - COS 0) - 

Km y '• J 

-|-n^|l — 3ne^”^ -|- nsik 0 -f 2e”^ — 2 cos^ 0 — cos0 -|- e"^^cos0 -|- 3ncos0e^”^(cos0 — 1) — ncos0 
+e'^^ cos^ 0(4 cos^ 0 — 5) I -|-n|3 cos^ 0 — 4 cos^ 0 + e“"^^(cos^ 0 — ^) -|- cos 0(1 + e”"^^) 


n(3)o = 0. 


(42) 
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For the solution given by Eq. (25) the basic vector has the non-vanishing components 


Ti = (^=ae“^ — \/2)(sin^ 0 -|- cos0), T 2 = (ae'^^ — \/2), 

v 2 v 2 


( 43 ) 


and the non-vanishing components of the axial vector held are 


= -n2V2e-^P + + 1){V2 - a)] {sin^ 0-t-cos 0} - 2e^(^-“)(\/2a - j . (44) 


A=-| 

3' 


[2^26-“^ + + 1)(^2 - a)] {sin^ 0 + cos 0} - (^20 - j e^. (45) 

The corresponding extra Hamiltonian [7] is given by 


6H = —ha = —- 


2 ^26-“^ + + 1)(V2 - a)] {sin2 0+cos0}-2e^(^-“H^a-e"“Oje.-c^, 

(46) 


and the component of the energy density has the form 


n(0)0 ^ ^ - ^) + COS 0(0^6“^ “ ^ - a) j . (47) 

The non-vanishing components of the momentum density have the form 
n{i)o = — a} 4- e““^(l — a\/2 — 3\/2 cos^ 0) -|- a(3 cos^ 0 — 1)) 

K, ^ ' 

pap / 

= — e““^(2\/2 cos 0(1 -1- cos 0) — \/2) -|- \/2e“^(8a^ — a cos^ 0 — 3a cos 0 -f 4a cos"^ 0 — 2a) 

2k, y 

-|-3a^e^“^(sin^ 0 + cos 0) + 2 cos^ 0(3 — 4 cos^ (j) — a? — 2a) + 2 cos 0(1 — a — a^) -|- 2a(a — 1)^ , 
= 0. (48) 


4. Main results and discussion 


We have applied the axisymmetric tetrad held given by Eq. (20) with three unknown 
functions of p and to the held equations (7). We have obtained two special solutions by 
taking particular values of the unknown functions E, F and K. The hrst one is given by 
Eq. (24) and is known as electromagnetic generalization of the Gddel solution [22], while the 
second one is given by Eq. (25) and is known as Gddel solution [22]. These solutions are 
special solutions of the non-linear P.D.E. given by Eq. (23) [27, 22]. 
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We have calculated the basic vector T“ dehned by Eq. (6) for the tetrad held (20). 
The components we have obtained as given by Eq. (27) are different from the components 
obtained in Ref. [27] which are 

T, = ^{F-e^)-k, n = -K'. (49) 

We have calculated the axial vector part applied Eq. (20) to Eq. (8). Our results are 
completely different from the results obtained in Ref. [27]. Therefore, the analysis of the 
extra Hamilton used in [27] is again discussed. The extra Hamiltonian is now recalculated 
for the tetrad field of Eq. (20) and is given by Eq. (33). 

The energy density and momentum density are calculated using the energy-momentum 
tensor of TEGR which is coordinate independent. For the energy momentum density and 
momentum density derived in Ref. [27] we have some comments: 

i) The energy momentum density and momentum density, are not correct. 

ii) In Ref. [27] a discussion is given for the case when the arbitrary parameter A = 1. It is 
not clear whey A = 1 is discussed in spite that for this case still the theory of Mpller deviates 
from GR. 

iii) When the tetrad (21) of Ref. [27] applied to the superpotential (12) of [27], it is logic to 
obtain the energy density and momentum density with the arbitrary parameter A. If we want 
to compare the result with GR we must take A = 0 not A = 1! Here in this work we have 
done the calculations of the energy density and momentum density using the gravitational 
energy-momentum tensor which is coordinate independent since it is constructed out from 
Hamilton structure. 
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